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We present a calculation of critical phenomena directly in continuous dimension d employing an 
exact renormalization group equation for the effective average action. For an Ising-type scalar field 
theory we calculate the critical exponents v{d) and 77(d) both from a lowest-order and a complete 
first-order derivative expansion of the effective average action. In particular, this can be used to 
study critical behavior as a function of dimensionality at fixed temperature. 
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INTRODUCTION 

One of the important quantitative successes of statistical 
physics is the precise understanding of critical phenom- 
ena based on the renormalization group Universal 
properties near second-order phase transitions are found 
to be insensitive to the details of the underlying the- 
ory and only to depend on symmetry, field content and 
the dimensionality of space d. The crucial dependence 
of critical behavior on d can be readily observed from 
the fact that no phase transition at nonzero tempera- 
ture or, equivalently, no corresponding renormalization 
group fixed point exists for d = 1. For higher dimensions 
phase transitions can occur, depending on the model. 
For instance, simple Ising-type one-component models 
can exhibit an infinite number of multi-critical points in 
d = 2 @ , a critical "Wilson-Fisher" or a tri-critical fixed 
point in d = 3, and a trivial "Gaussian" fixed point for 
d = A. 

Continuity in the dimensionality of space is a cru- 
cial ingredient for some of the standard methods for 
the description of critical phenomena. In particular the 
extrapolation from d = 4 to lower dimension (4 — 6 
expansion Q) is one of the most frequently applied meth- 
ods for studying the non-perturbative dynamics near a 
non-trivial fixed point of the renormalization group equa- 
tion. In this language the strong infrared divergences en- 
countered by perturbation theory at fixed dimension are 
generated by an expansion around the "wrong" (Gaus- 
sian) fixed point. The e-expansion is used to interchange 
limits which allows one to follow "perturbatively" the 
relevant fixed point. The involved extrapolation in the 
dimension d is a crucial ingredient for this procedure. In 
particular, the typical expansion parameter e = 1 is not 
small and improvement procedures such as Borel trans- 
formation and conformal mapping are required to obtain 
accurate results 3J. A direct calculation in continuous 
dimensions would provide important insight concerning 
the validity of such an approach. 

In this work we calculate critical phenomena for con- 
tinuous dimension d within a one-component scalar field 



theory with quartic self-interaction. Computations in ar- 
bitrary dimension can be performed 0, using the effec- 
tive average action [(J, |7j , which is infrared regulated 
by a momentum cutoff. The change of Tk with the scale 
k is described in terms of an exact flow equation 0|: 
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where T k '[<p] is the exact inverse propagator and Rk de- 
notes the cutoff function. We investigate the theory in 
complete first-order in a derivative expansion, where 



r, = d 



(2) 



with p = qb 2 /2. Critical phenomena can be described 
from the flow of the most general effective potential Uk(p) 
consistent with the symmetries, as well as the wave func- 
tion renormalization Zk (p) Q • This approach allows one 
to perform direct calculations for a given continuous di- 
mension parameter d. In particular, this can be used to 
study critical behavior as a function of d at fixed tem- 
perature. 



EVOLUTION EQUATIONS 

For the truncation J2J) the exact renormalization group 
equation Q is reduced to partial differential equations 
for U}Ap] and Zk(p)- These equations can be found in 
Refs. pUsf- For the cutoff function we employ |S( 



R k = z k (k 2 -p 2 )e(k 2 - P 2 ) 



(3) 



with Zk = Zk(p = po) evaluated at vanishing external 
momentum. The flow of the function Z k is determined 
by the anomalous dimension rj = — 9 t lnZ^. Our quanti- 
tative results are based on a complete numerical solution 
of the equations for the first-order derivative expansion. 
For comparison, we solve in addition a simpler approxi- 
mation with a field independent wave function renormal- 
ization Zk- For our analytical discussion we concentrate 
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on the latter, for which the flow of the rescaled effective 
potential u(p) = k~ d Uk as a function of the renormalized 
and dimensionless field variable p — Z k k 2 ~ d p is described 
by 



dtu = —du + (d — 2 + rj)pu 
4«d ( , T) 
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with vj 1 = 2 d+1 n d ^T{d/2) and u' = du/dp. Here 
t = ln(A;/A), where A denotes a high momentum scale. 
Most of the relevant qualitative properties can be ob- 
tained from the /c-dependence of the location k of the 
minimum of u(p). By taking the total t-derivative of 
u'(p = k) = the scale-dependence of the minimum for 
k > is found to be: 



d t K = P K = -(d-2 + r])K + 
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Nr = ii-^v+ 2K ?( ( y. (7) 



d+2J (1 + 2kA) 2 ' 

Here the couplings are related to p-derivatives of the 
rescaled potential evaluated at the minimum, i.e. A = 
u"{k) and u 3 — u"'(k). Note that (3 K is parametrized 
in terms of the couplings A and u$. Implicitly, however, 
/3 K depends on an infinite number of couplings since the 
correponding /3-function for A depends on the higher cou- 



plings U3 and 114 



(k), the one for U3 depends on U/± 



and U5 = u^ 5 >(k) etc. 

A second-order phase transition is characterized by a 
fixed point or scaling solution for which d t u' = and, 
in particular, (3 K — 0. All renormalized dimensionless 
couplings become independent of the scale k and take on 
their "fixed point values" «*, A*, M3* etc. Typically, n 
corresponds to the relevant parameter such that critical 
behavior can be observed only for a fine-tuned initial 
value k(A) ~ k c with (k(A) — k c ) ~ (T c — T). If k is the 
only relevant parameter the flow of all other couplings 
A, U3 and higher derivatives of the potential follow fixed 
"critical trajectories" A(k, d), u 3 (n, d) — so-called partial 
fixed points. (More precisely, this holds if T is in the 
vicinity of the critical temperature T c and k is sufficiently 
below A.) Inserting the critical trajectories, j3 K becomes 
only a function of n and d: 



P k (k, d) = (3 k (k, X(k, d),...,d). 



(8) 



For arbitrary d the critical behavior is then directly linked 
to the properties of (3 K (K,d) — in particular its zeros. 
The precise form of the critical trajectories X(k, d) etc. is 
often not known analytically but it can be determined nu- 
merically. However, we will see that many crucial prop- 
erties of (3 K (k, d) do actually not involve a very precise 



analytical understanding of the critical trajectories. We 
emphasize that the use of A(k, d) etc. in the flow equa- 
tions is valid only in the vicinity of the critical trajectory. 
The analytical discussion in the next section is therefore 
limited to the behavior close to criticality. The full nu- 
merical treatment below is not restricted to this. 



ASYMPTOTIC BEHAVIOR FOR LARGE k 

As a general feature, (3 k {k, A, d) is analytic in all 
arguments except for A = 0. This does, however, not 
imply automatically that /3 k (k, d) is analytic since the 
solutions for the critical trajectories of X(k, d) and the 
higher couplings can be non-analytic. For fixed K, A, etc. 
the function /3 K decreases with increasing d. According 
to and © one has 

12v d 



a. 

To discuss (3 K for n —. 
havior is given by 

r (2 - d)n 
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for k = . 
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00 we note that the leading be- 
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00, d > 2 
00, d < 2 
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As is shown below, here Coo is 3- constant which de- 
pends on d. In particular, in two dimensions one finds 
clo* = — l/(47r) ||. For positive c£? the asymptotic be- 
havior can be related to a so-called "zero-temperature 
phase transition" . In this case critical behavior can be 
associated with a fixed point for dtK" 1 = —k~ 2 /3 k at 
1/k = 0. 

To establish the leading behavior (|10|l for d<2we first 
note that the region of large n admits slowly evolving 
solutions if 77 is near 2 — d. From Eqs. © and © one 
observes 
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d(d + 2)^ \d-2 + ip 

which shows that the presence of fixed points in the re- 
gion of large k 3 A 2 requires values of 77 close to 2 — d. We 
parametrize 77 = 2 — d + C/ K > where £ is a K-independent 



constant. Introducing the notation u(p) 
finds from Eq. (0J: 

dth Is 
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Since we are interested in values for p near k we consider 
the variable x = p — k and neglect terms suppressed by 
powers of x/k as well as £/ ' k, i.e. 



d t h \ x 



d t h \p 
-dh 4 



hh'd t K 
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We observe that k does not appear anymore in ltl^|) 
and consequently quantities like h 2 = h"(x — 0), h 3 = 
h"'(x — 0) are not expected to depend on n. In leading 



order for 



oc we therefore have A 



-1/2 



k^ 1 / 2 and Ak ~ k 1 / 2 » 1 as well as nu 3 ~ k 1 / 2 » A. In 



particular, N r is also independent of k and given by 

d h 3 



N r = 



3(d + 2] 



h\ 



This shows that the leading behavior of p K for n 
and d < 2 is described by 



(14) 



(15) 



Further information can be obtained from the scaling 
solution dth\ x = 0. (Note that the scaling solution for 
h does not correspond to an exact fixed point since n 
evolves according to i|15[l.l At the potential minimum we 
have h'(x = 0) = up to corrections ~ 1/k. From ijT3*|) 
one then finds 

112 = d(d%h ■ (16) 

The integration constant h Q = h(x = 0) has to be chosen 
such that the lowest order relation for the anomalous 
dimension, ijq = 2 — d, holds. According to Eq. (J3J) at 
this order 

18v d: 
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which fixes N r as a function of h 2 : 



N r 



(2 - d)d 



2v d 3(d + 2)/i 2 



(17) 



(18) 



and h 3 by i|14[) . Finally, h 2 is constrained by the require- 
ment that the scaling solution should extend to the whole 
range —k<x<k. For d — > 2 one concludes from l)18|) 
that N r — ► 0. As a consequence, h is determined by the 
motion of a particle in a logarithmic potential, i.e. for 
d — > 2 the lowest order scaling solution for l|13l) reads 



{h'f 
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(19) 



d(d + 2) 

It is striking that the simple equations Q , J5| and © 
reproduce correctly the low-temperature behavior of the 
exact Ising model for d = 1, which is characterized by 
7] = 1 and essential scaling: For T — > from the leading 
behavior i|15fl one has 



(20) 



Associating the scale k s where n(k s ) = to the inverse 
correlation length up to a factor one finds 



£A ~ exp 



(^ 



) ) 



exp(2J/T) 



(21) 



with k(A) = b/T, b ■■ 



2 Cqo J • 



NON-ZERO TEMPERATURE PHASE 
TRANSITION 

From the asymptotic behavior discussed above one can 
draw a number important conclusions about the presence 
of a phase transition at non-zero temperature, character- 
ized by a finite value . We note that a positive value of 
(3 K for k = and negative values for large k require the 
/3-function to vanish at an intermediate k — k*. Since 
clo* = — l/(47r) Eq. (|10|l immediately tells us that for all 
d > 2 a finite value exists for which j3 K = 0. For 
instance, in d = 3 the corresponding phase transition is 
described by the well-known Wilson-Fisher fixed point 
in the presence of only one relevant parameter 

From the exact solution of the one-dimensional Ising- 
model one has Co$ > and (3 K is above a positive constant 
for all k. Therefore, no fixed point with a finite oc- 
curs in d = 1. If the /3-function depends continuously on 
d we conclude that the zero-temperature transition and 
the absence of a fixed point with finite extend to a 
finite range of d > 1. (It is crucial in this respect that 
the limiting value is strictly positive.) As a conse- 
quence, the lower critical dimension, below which phase 
transitions at nonzero temperature are absent, has to be 
larger than one. 

The lower bound on the critical dimension, d\ ow > 1, 
can be supplemented by an upper bound using similar 
continuity arguments. For this we note that the cor- 
relation length exponent v at an infrared-unstable fixed 
point is given by the first derivative of /3 k (k, d) at =/= 0: 



d(3 K (n,d) 
dn 



(22) 



For finite v this yields the usualpower law for the cor- 
relation length, £ ~ \T — T c \~ v [3j. For a vanishing first 
derivative the exponent v diverges. In this case the power 
law is replaced by essential scaling. A finite value for v 
corresponds to nonzero negative df3 K /dn\ K ^ < 0. Since 
the derivative of the /3-function (and therefore v) de- 
pends continuously on d, we conclude that for finite v 
at a given d also a finite neighborhood in d exhibits a 
phase transition. At the lower critical dimension, below 
which no phase transition with a finite occurs, either 
— > oo or v must diverge and d\ ow is characterized 
by essential scaling. The finitcness of v for d — 2 follows 
from universality and the exact value for the Ising model, 
v = 1 3- Therefore v must also remain finite in a finite 
range d < 2. This continuity argument therefore implies 

rflow < 2. 

We emphasize that for a second-order phase transi- 
tion, apart from a vanishing /3 K -function, all couplings 
have to take on fixed-point values and the critical tra- 
jectories X(k, d) etc. have to be well defined. In partic- 
ular, the character of a fixed point is influenced by the 
stability properties with respect to the other couplings. 
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FIG. 1: The critical correlation length exponent v and anoma- 
lous dimension rj as a function of the dimension d for two dif- 
ferent approximations: complete first-order derivative expan- 
sion (open circles) and for a field-independent wave function 
renormalization (full circles). We also show the results of the 
e-expansion up to e 5 (solid line) [8], which are plotted for the 
range < e < 3 for illustration, along with available exact 
results (open squares). 



To establish this we turn to the numerical solution of 
the full first-order derivative expansion, as well as for 
the approximation with a field-independent wave func- 
tion renormalization as discussed above. Since the ap- 
proximations include a most general potential term u(p) 
the description takes into account an "infinite number of 
couplings" u n = u (n) (p)|p =K . 

We numerically solve the evolution equations, which 
are given e.g. in Eqs. (3.8) and (3.10) of Ref. [llj, for the 
one-component theory to first order in the derivative ex- 
pansion. The equations were discretized on a grid for the 
variable p with varying grid-size (typically thirty points 
and much higher resolution in lower dimensions). The re- 
sulting ordinary differential equations were solved using 
a fifth-order Runge-Kutta algorithm and a fourth-order 
Cash-Karp stepping routine. The relative local errors of 
the integration routine were below 10 -6 . Results for sam- 
ple exponents in various dimensions were cross-checked 
with different numerics (in particular, the integration of 
the flow equation of the potential u(p) instead of its first 
derivative u'{p), as well as a fixed grid instead of run- 
ning discretization points ~ k). The quantitative differ- 
ences in the exponents using the different numerics were 
manifest only in the third to fifth significant digit. 

The results for the critical exponent v(d) and the 



anomalous dimension 77(d) are plotted in Fig. ^ for the 
two different employed approximations. For comparison 
we also include the results obtained from the expansion 
around 4 — e dimensions up to e 5 0, as well as avail- 
able exact results. The results from the e-expansion have 
been plotted for values as large as e — 3 for illustrational 
purposes though it is typically employed only for e = 1. 
In contrast, we emphasize that our results correspond to 
direct calculations for the corresponding dimensions. In 
particular, they do not involve any improvement proce- 
dure such as Borel transformation and conformal map- 
ping employed for the expansion around 4 — e 0, . 

To numerically obtain v(d) we calculate a series of 
near-critical trajectories with varying (/c(A) — k c ) ~ 
(T c — T). These calculations become time intensive for 
low dimensions, and grow substantially for the more so- 
phisticated approximation with a field-dependent (p) . 
This is why we have computed v(d) only down to about 
d = 1.5 in Fig. ^ In contrast, the anomalous dimension 
77 can be directly infered from the critical trajectory and 
is presented in Fig. ^ for even lower dimensions. 

As d is lowered from the upper critical dimension d — 4 
both v{d) and 77(d) grow. One observes a very good agree- 
ment of the results from our first-order derivative expan- 
sion and the e-expansion results down to about d = 3. At 
d = 2 we can compare with the exact values u = 1 and 
77 = 0.25. The renormalization group results for v with 
the O-cutoff function Q show a rather strong depen- 
dence on the approximation around d — 2. The strong 
increase of v{d) as d is further lowered sets in much ear- 
lier for the approximation with a field-independent wave- 
function renormalization. Interestingly, for the anoma- 
lous dimension the dependence on the approximation is 
diminished for lower dimensions. In particular, we ob- 
serve that already the less sophisticated approximation 
approaches the exact one-dimensional result. This may 
be related to fact that in one dimension the propagator 
shows a simple quadratic dependence on momenta for 
77= 1 [3. 

CRITICAL BEHAVIOR AS A FUNCTION OF 
DIMENSIONALITY 

We point out that critical behavior can also be studied 
as a function of d for fixed k(A) or, equivalently, temper- 
ature T. In the critical region we find that for fixed T 
the correlation length varies with d as 

t~\d c -d\- v . (23) 

Here d c is the dimension for which T becomes the crit- 
ical temperature and v = v(d c ) is the same critical ex- 
ponent as computed above for the variation with T for 
fixed d. This property can be understood as a simple con- 
sequence of continuity. Assume some fixed point where 
/3 K (K*{d c ),d c ) = and (d/3 K /dK)(K*(d c ),d c ) = -v~ l (d c ) 
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with finite critical exponent v(d c ). Starting with criti- 
cal initial values («*(<i c ),d c ) at some scale A the correla- 
tion length £ diverges. For initial values (k„ (d c ) — A, d c ) 
we have the well-known critical behavior (A > in the 
symmetric phase) £(A,d c ) ~ A~"( dc K Continuity tells 
us that for d = d c — A c i there is a new fixed point 
with /3 K (fc*(d c — A),d c — A) = 0. Instead of varying 
A at fixed d c one may keep a fixed k(A) = K*(d c ) and 
vary A^. With n^(d c ) = K*(d c — A^) — A this yields 
£(0,4 - A d ) ~ A- u( - d "- Ad l (Note that A d > for the 
symmetric phase.) For A^ — » we can relate A to A d 

A = -^(d c )A d = -y(d c )^{K*{d c ), d c ) A d . (24) 

Generically one has d(3 K /dd > and we conclude that the 
dependence of £ on d at the critical point (d c ) is char- 
acterized by the standard critical exponent v = v(d c ), 
which agrees with the results from the numerical solu- 
tion to very good accuracy. This is investigated for the 
corresponding A-component model in Ref. |l2j . We re- 
call, however, that this argument is only valid for finite 
and away from essential scaling at d\ ow . 



CONCLUSIONS 

We have presented a direct calculation of critical phe- 
nomena in continuous dimension emloying a complete 
first-order derivative expansion. The results for the crit- 
ical exponents show good agreement with those found 
from the e-expansion for not too small dimensions down 
to about d = 3. The deviations grow substantially for 
lower dimensions. To get an idea about the reliability of 
our results we have re-calculated the exponents using a 
less sophisticated approximation with a field-independent 
wave-function renormalization. It is striking to observe 
that though the critical behavior depends comparably 
strongly on the employed approximation scheme for in- 
termediate values of (i, this has only a small effect on 
some nontrivial properties at sufficiently low dimensions: 
even the less sophisticated approximation is sufficient to 
approach the exact result for the anomalous dimension 
rj = 1 in d = 1. It has recently been demonstrated for the 
three-dimensional one-component model that the next 
order (<9 4 ) of the derivative expansion leads to signifi- 
cantly improved results [l^j ] compared with lower orders 
calculations. It would be very interesting to investigate 
the apparent convergence of the derivative expansion for 
lower dimensions as well. 



We have pointed out that with the employed methods 
one can use the dimensionality d as a "relevant" parame- 
ter to study deviations from criticality. For fixed temper- 
ature T the deviations from the critical dimension d c are 
described by the same universal critical exponents which 
parametrize the deviations from T c for fixed d. So far we 
are not aware of any clear experimental application of the 
calculated critical phenomena in fractional dimensions. 
However, an effectively similar behavior may be observed 
e.g. for the magnetic properties of thin films, showing 
a dimensional crossover from three- to two-dimensional 
scaling properties by varying the film thickness [T^j . 
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